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FOREWORD 

This  report  documents  an  evaluation  of  statistical  fracture  theories 
as  applied  to  aerodynamic  heating  of  full-scale  seeker  domes.  Dome 
fracture  data  obtained  from  aerodynamic  heating  tests  of  magnesium 
fluoride  domes  were  used  to  evaluate  these  fracture  theories.  The  work 
was  performed  during  fiscal  years  1978  and  1979. 
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INTRODUCTION 


The  launch  envelopes  of  heat-seeking  missiles  used  on  current  combat 
aircraft  are  frequently  restricted  by  the  thermos true tural  limits  of  the 
infrared  (IR)  dome  used  in  the  seeker  system.  These  limits  are  derived 
from  the  predicted  fracture  of  the  ceramic  materials  used  in  IR  dome 
construction.  In  general,  ceramics  are  brittle  materials  that  exhibit  a 
wide  range  of  fracture  strengths  under  apparently  identical  loading  con¬ 
ditions.  Consequently,  unrealistically  severe  missile  launch  envelope 
restrictions  are  imposed  on  the  aircraft  in  order  to  compensate  for  the 
uncertainties  in  the  dome  fracture  strength. 

The  results  of  a  previous  investigation^  gave  indications  that  better 
definition  of  such  launch  envelope  restrictions  could  be  achieved  through 
the  use  of  the  statistical  fracture  theories  of  Weibull^jS  and  Batdorf. ^*5 
The  limited  amount  of  IR  dome  fracture  data  extant  at  the  time  of  this 
previous  investigation  prevented  a  comprehensive  evaluation  of  these 
fracture  theories.  In  addition,  the  mathematical  details  required  to 
utilize  material  fracture  statistics  from  a  variety  of  fracture  strength 
test  methods  were  not  yet  developed  for  the  Batdorf  theory.  The  current 
investigation  extends  the  usefulness  of  the  Batdorf  theory  bv  developing 
these  mathematical  details  for  three  commonly  employed  strength  test 
methods . 


1  W.  R.  Compton.  Application  of  Statistical  Fracture  Criteria  to  the  Problem  of  Predicting  Infrared 
Dome  Thermal  Shock  Failures.  China  Lake.  Calif..  NWC.  January  1978.  (NWC  TP  6010.  publication 
UNCLASSIFIED.) 

-  w.  Weibull.  “Statistical  Theory  of  Strength  of  Materials,”  Ing.  Vetenskaps  Akad.  Hand!..  No.  151. 
45  pp  (1939);  Ceram.  Abstr. ,  19  [3]  78  (1949) 

3  w.  Weibull.  “Statistical  Distribution  Function  of  Wide  Applicability,”  /  Appl.  Meek.  1S[3]293 
(September  1971) 

4  The  Aerospace  Corporation.  Fracture  Statistics  of  Brittle  Materials  with  Intergranular  Cracks .  by 
s.  B.  Batdorf.  El  Segundo,  Calif.,  10  October  1974.  (SAMSO-TR-74-210.  publication  UNCLASSIFIED.) 

5  - .  Fracture  Statistics  of  Isotropic  Brittle  Materials  with  Surface  Flaws,  by  S.  B.  Batdorf.  El 

Segundo,  Calif.,  3  December  1973.  (SAMSO-TR-73-378.  publication  UNCLASSIFIED.) 
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Fracture  statistics  for  IR  dome  materials  are  usually  determined 
from  flexural  tests  of  small  laboratory  material  specimens.  Such  material 
specimens  frequently  do  not  have  the  same  internal  grain  structure  or 
external  surface  preparation  that  exists  for  an  IR  dome.  Consequently, 
any  difference  in  the  method  of  fabrication  between  the  small  material 
specimens  and  an  IR  dome  may  yield  erroneous  fracture  statistics  for  the 
dome.  Two  sets  of  fracture  statistics  for  the  material  magnesium 
fluoride  (MgF2)  were  used  in  this  investigation  in  order  to  explore  this 
problem  area.  One  set  of  fracture  statistics  was  obtained  from  four- 
point  bending  tests^  of  specimens  with  surface  scratches  less  than 
6(10)"^  inch  deep  (15  pm)  and  another  was  obtained  from  concentric  ring 
bending  tests?  of  specimens  with  surface  scratches  less  than  4(10)*^  inch 
deep  (1  pm).  Both  sets  of  material  test  specimens  were  fabricated  from 
MgF2  plate  stock. 

The  objective  of  the  present  investigation  was  to  obtain  a  comprehen¬ 
sive  evaluation  of  the  aforementioned  fracture  theories  by  direct  com¬ 
parisons  with  full-scale  IR  dome  fracture  statistics.  The  fracture 
statistics  for  IR  domes  were  obtained  by  subjecting  thirty  MgF2  domes  to 
the  same  simulated  free- flight  thermal  environment  of  the  NWC  T-Range 
hot  gas  facility.  The  statistical  fracture  theories  were  evaluated  by 
comparing  the  predicted  and  observed  probability  of  dome  fracture  as  a 
function  of  time  after  the  simulated  missile  launch. 


DOME  FRACTURE  TESTS 


Fracture  data  for  IR  domes  were  obtained  at  the  NVC  T-Range  hot  gas 
facility  by  subjecting  30  MgF2  domes  (Figure  1)  to  a  thermal  environment 
simulating  the  free- flight  profile  shown  in  Figure  2.  The  T-Range 
facility  utilizes  an  axial  flow  propane  burner  in  a  blowdown  wind  tunnel 
to  produce  free-f light  aerodynamic  heating  levels.  The  MgF2  domes  were 
mounted  in  the  free-jet  issuing  from  a  convergent-divergent  nozzle  as 
shown  in  Figure  3.  The  free-flight  thermal  environment  was  simulated  by 
varying  the  total  pressure  and  temperature  histories  in  the  T-Range 
facility  in  such  a  way  as  to  match  the  total  conditions  on  the  dome  during 
tree-flight  (Table  1).  The  range  of  total  pressure  and  temperature  pro¬ 
files  actually  observed  during  the  test  series  is  shown  in  Figures  4  and  5. 


^  Naval  Weapons  Center.  Report  of  Modulus  of  Rupture  Tests  on  Magnesium  Fluoride  IR  Dome 
Material,  by  R.  L.  Smith.  China  Lake,  Calif.,  NWC.  8  August  1972.  (Reg.  Memo  4062-010-73. 
UNCLASSIFIED.) 

7 - .  Unpublished  test  data  for  fracture  strength  of  magnesium  fluoride  discs  subjected  to 

con;:ntnc  ring  ioadinc.  by  R.  J.  Schiltz  and  G.  A.  Hayes.  China  Lake.  Calif..  NWC.  l0"5. 
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IR  DOME  CALORIMETER  DOME 


MAGNESIUM  FLUOR  IOE  304  ST  /UNLESS  STEEL 


FIGURE  1.  Magnesium  Fluoride  Dome  Geometry. 


The  variation  of  the  aerodynamic  heat  transfer  coefficient  with  time 
and  dome  streamwise  position  was  determined  from  the  thermal  response  of 
a  calorimeter  dome  constructed  from  0.060  inch  thick:  304  stainless  steel 
(see  footnote  1).  The  heat  transfer  coefficient  variation  with  time  and 
dome  streamwise  position  is  shown  in  Figures  6  and  7.  The  heat  transfer 
levels  obtained  in  the  free  jet  were  found  to  be  more  severe  than  those 
encountered  during  free  flight.  This  was  due  to  the  high  turbulence 
levels  inherent  in  the  T-Range  facility.  The  stagnation  point  heat 
transfer  coefficient  was  roughly  40*  higher  than  expected  during  a 
missile  free  flight.  In  addition,  boundarv  layer  transition  occurred 
much  closer  to  the  dome  stagnation  point  than  would  be  expected  during 
free  flight. 
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NOTES:  ALTITUDE  -  30,000  FT  (CONSTANT) 
STANDARD  DAY  ATMOSPHERE 


TIME  AFTER  LAUNCH  ,  SEC 


FIGURE  2.  Free-Flight  Profile. 
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FIGURE  3.  IR  Dome  Test  Setup. 


TABLE  1.  Simulated  Free-Fl ight  Conditions. 


Free-flight  conditions 

T-Range 

nozzle 

exit  conditions 

c, 

Mach 

— 

P  , 

T  , 

V 

T 

D 

k  r  * 

sec 

no . 

t 

psia 

c 

°R 

psia 

3R 

psia 

0 

2.25 

51 

S28 

31 

2.54 

64 

828 

31 

1 

2.64 

93 

981 

41 

2. 54 

86 

981 

41 

2 

3.06 

178 

1170  1 

55  ! 

2.53 

114 

1170 

55 

3 

3.49  1 

336 

1386 

71  i 

2.  52 

147 

1386 

71 

4 

3.90 

602 

1613 

88 

2.51 

1S2 

1613 

88 

5 

4.20 

910 

1791 

102 

2.50 

210 

1791 

102 

6 

3.84 

■ 

554 

1579 

85 

i 

2.51  1 

176 

1579 

35 

-  Total  pressure  upstream  of  normal  shock 

*  *  Total  pressure  downstream  of  normal  shock 

^  n 

u 

T  =  Total  temperature 
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6.  DEGREES 

FIGURE  7.  Streanvise  Heac  Transfer  Coefficient  Distribution 


Temperature  and  thermal  stress  levels  within  the  M<jF2  dome*  during 
the  simulated  free  flight  were  predicted  from  thermal  and  stress  computer 
models®*9  (see  also  footnote  1),  and  the  observed  total  temperature  and 
heat  transfer  coefficient  profiles*  The  accuracy  of  the  computer  models 
have  been  verified  experimentally  in  an  earlier  investigation  (footnote  1) 
The  predicted  thermal  stress  profile  corresponding  to  the  time  variation 
in  facility  total  pressure  and  temperature  is  shown  in  Figure  3  for  the 
dome  stagnation  point.  The  predicted  strearawise  variation  of  thermal 
stresses  in  the  dome  at  the  end  of  the  free-flight  boost  phase  is  shown 
in  Figure  9.  It  was  found  (see  footnote  1)  that  the  metal  dome  housing 
could  be  idealized  by  a  simply  supported  dome  rim  with  small  less  of 
accuracy.  The  predominant  effect  of  the  housing  on  the  dome  stress  dis¬ 
tribution  was  to  change  the  dome  temperature  distribution  locallv. 

All  of  the  MgF2  domes  fractured  during  the  free-flight  aerodynamic 
heating  stimulation.  Time  to  dome  fracture  after  *-he  simulated  missile 
launch  was  obtained  from  400  frame/second  camera  records  of  each  dome 
test.  The  probability  of  MgF?  dome  fracture  as  a  function  of  time  after 
simulated  missile  launch  is  shown  in  Figure  10.  All  domes  were  subje'red 
to  a  thermal  environment  corresponding  to  the  Mach  2.25  launch  condition 
for  a  minimum  of  20  seconds  prior  to  the  start  of  the  simulated  free  fiign 
Mote  that  50%  of  the  domes  failed  within  4.5  seconds  (Figure  10)  at  a 
corresponding  stress  of  less  than  10,000  psi  (Figure  8).  The  published 
value  of  the  mean  fracture  stress  obtained  from  small  material  specimens 
is  on  the  order  of  20,000  psi.^ 

WEI BULL  FRACTURE  THEORY 

The  fracture  theorv  developed  by  Weibull  (see  footnote  21  assumes 
that  fracture  occurs  when  the  tensile  stress  exceeds  the  srrengtn  the 
weakest  flaw  present  in  a  material  sample.  It  is  further  assumed  ■  \a r  tho 
flaws  (type  unspecified)  occur  in  a  manner  such  that  the  nrobab i i i i  *•  of 
fracture  for  the  brittle  material  is  described  by  the  expression: 


®  Naval  Weapons  Center.  Aerodynamic  Heating  of  Spherically  Tipped  Cylinders ,  Cones  and  Ogives 
Using  the  General  Thermal  Analyzer  SIX  DA ,  by  W.  R.  Compton.  Cluna  Lake.  Calif..  WVC.  June  !9"4 
(\WC  TN  4061  *172.  publication  UNCLASSIFIED.) 

^  The  Aerospace  Corporation.  SAAS-III  Finite  Element  Stress  Analysis  of  A  Asymmetric  ano  Plane 
Solids  with  Different  Ortho  tropic.  Temperature-Dependent  Material  Properties  on  Tension  ami  C  impress- 
sion .  by  R.  M.  Jones  and  J.  G.  Crose.  Los  Angeles.  Calif.,  June  1971.  (TR-0059  (56810-531-1.  publication 
UNCLASSIFIED.) 

10  The  Eastman  Kodak  Company.  Kodak  Irxran.  Infrared  Optical  Materials.  197].  (Kodak  Publication 
U*72,  UNCLASSIFIED.) 
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FIGURE  10.  Experimental  Probability  of  Dome  Fracture. 
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where 


=  probability  of  fracture 
0  =  applied  tensile  stress 


a  =  tensile  stress  at  which  P-  *  0 
u  f 

0q  =  a  scaling  parameter 
m  =  a  flaw  density  parameter 


H  *  either  the  surface  area  or  the 
volume  of  the  material 


The  parameters  0U,  a0,  and  m  are  determined  by  a  least-squares  fit  of 
the  fracture  strength  data  obtained  from  flexural  strength  tests  (see 
footnote  3) .  The  values  of  0Q  and  m  for  MgF2  have  been  determined 
using  Weibull’s  two-parameter  (c  =  0)  method  using  the  expressions  below. 


N  N  N 

V  X  T)  Y  -  M  52  X  Y 

"  n  "  n  "  n  n 
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n=l 
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n-l 


2  >J 

- £  * 

n=l 


NWC  TP  6226 


where 


H 

L 

V 

w 

X 

n 

Y 

n 


material  specimen  thickness  parallel  to  the  applied  load 
(Figure  11) 


material  specimen  length  between  supports 
material  specimen  volume  between  supports 
material  specimen  width  normal  to  the  applied  load 
Inc 


n 


£n£n 


n  +  1 
n  +  1  -  n 


n  =  number  of  material  specimens  that  failed  at  a  maximum  tensile 
strength  of  or  lower 

N  *  total  number  of  material  specimens  tested 

a  «  maximum  tensile  stress  at  fracture  for  the  n11*1  material  specimen 
r 


probability  of  survival 


UNIAXIAL  TEST  SPECIMEN 


LOAD 

t  * 

□z 


UNIFORM 
RING  LOAD 


EQUIB1AXIAL  TEST  SPECIMEN 

FIGURE  11.  Material  Specimen  Geometry. 
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A  summary  of  values  for  the  Weibull  parameters  obtained  for  MgF  7  is 
given  in  Table  2.  The  fracture  data6  are  the  data  obtained  from  flexural 
strength  tests  of  material  specimens  with  surface  scratches  less  than 
6(10)"^  inch  (15  pm)  deep.  A  comparison  of  the  results  of  the  least- 
squares  approximation  and  the  fracture  strength  data  is  shown  in 
Figures  12  through  16. 

1  The  probability  of  fracture  for  a  given  dome  geometry  is  determined 

as  follows: 

1.  The  temperature  distribution  within  the  dome  is  determined  as  a 
function  of  time. 

i 

k 

2.  The  thermal  stresses  within  the  dome  are  determined  at  a  partic¬ 
ular  time  from  the  temperature  distribution  and  a  finite-element  stress 
analysis . 

V 

«,  3.  The  probability  of  survival,  Ps  =  1  -  Pf,  is  determined  for 

<  each  volume  (surface)  element  (An)  in  the  stress  model  bv  using  the 

*  maximum  principal  stress  (tensile)  in  that  element  (a) ,  and  the  approxi- 

•  mat ion, 


4.  The  probability  of  dome  survival  is  then  determined  for  that 
particular  time  by  taking  the  product  of  the  individual  volume  (surface) 
element  probabilities  of  survival. 

Weibull1 s  method  is  widely  used  to  predict  brittle  material  fracture. 
There  is,  however,  no  theoretical  way  of  emonstrating  that  the  method 
is  valid  for  materials  that  exhibit  large  fracture  strength  dispersions. 
3atdorf  (see  footnote  4)  has  demonstrated  analytically  that  the  physical 
model  implied  by  Weibull 's  method  is  not  satisfactory  when  the  flaw 
density  parameter,  m,  is  greater  than  six. 


TABLE  2.  Uniaxial  MgF0  Fracture  Data. 


Temperature , 

°  F 

_  _  -  —  i 

Volume  theory, 

V  psi 

i 

Surface  theorv,  | 

t  m 

V  psl  ; 

_ i _  . 

75  1 

7,507 

19,814 

!  8.  5 

250 

10,152 

19, S97 

12.  73 

500 

8,625 

20,302 

9.77 

750 

4,470 

20,371  i 

5.  12 

1000 

2,279 

17,269 

_ j 

3.75 

probability  of  failure.  % 
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BATDORF  FRACTURE  THEORY 


VOLUME-DISTRIBUTED  CRACKS 

The  fracture  theory  developed  by  Batdorf  (footnote  4)  assumes  that 
the  flaws  existing  within  a  brittle  material  are  penny-shaped  intergranular 
cracks.  In  addition,  Batdorf  devised  a  means  of  accounting  for  the  effects 
of  crack  orientation  on  the  probability  of  fracture  in  triaxial  stress 
states.  Batdorf  arrived  at  the  following  expression  for  the  probability 
of  survival  for  a  small  volume  element  at  a  uniform  stress  state: 


where 


M  =  probability  that  a  crack  of  critical  strength  for  less)  is 
present  in  the  volume  element  (IV) 

—  =  probao ili ^y  that  the  crack  is  oriented  in  a  direction  such  that 
the  seres  normal  to  the  crack  plane  exceeds  the  critical 
crack  strength 

=  crack  critical  strength 


The  probability  that  the  crack  is  oriented  in  a  direction  that  will 
result  in  fracture  is  determined  from  the  expression, 


T 

"T 


cos6  d4> 


where 


ry  k  K  -  (K  -  K  )sin  : 

2  n  _  C  V  V  2  ,  ,  .  .  .  .  , 

cos  U  *  - — ^ -  (shear  insensitive  cracks) 

1-  K  -  (K  -  K  )sin*~  2 
y  v  z 
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FIGURE  17.  Volume-Distributed  Crack  Coordinate  System. 
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The  preceding  integral  has  been  evaluated  for  the  following  stress 
states: 

Shear  Insensitive  Cracks 


(K  »  0,  K  =0) 
v  z 


(K  *  1,  K  -  0) 

y  z 


(uniaxial ) 


(equib iaxial) 


*See  Appendix  A. 
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Shear  Sensitive  Cracks 


(K  •  0,  K  =0) 

y  z 

*  1  -  Kq  (uniaxial) 

(K  =  1,  K  =0) 

y  z 

=  ^1  -  (equibiaxial ) 

(K  »  1,  K  =1) 

V  2 

=  1  (equitriaxial ) 

Evaluation  of  the  integral  for  0.1^1  for  shear  sensitive  cracks  in  a 
triaxial  stress  state  involves  rather  messy  hypereliiptic  integrals  and 
is  best  treated  via  numerical  integration  techniques. 

If  the  critical  crack  strength  Parameter,  M,  is  expressed  as  a 
ser  ies , 


n=i  n=l 


then  the  probability  of  survival  can  be  determined  analytically  for  the 
following  special  cases: 


4nP 


-  nP 


-vV 

La  (2n  + 


Shear  Insensitive  Cracks 
_n 

(uniaxial,  pure  tension) 


n»l 


1) 


V  X  '  ^  ^  ^ 

"  2  2a  MTmIT+T)  (uniaxial,  pure  bending) 


n=l 


In? 


-  t  ^  ^  f)  B(n’f)  (equibiaxial,  pure  bending) 


n=l 
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where 


l r 

Shear  Sensitive  Cracks 
n 


?nP  *  -V 
s 


in? 


in? 


En  m 
(n  + 

n*l 

..  ■  a  a 

VV-_ILJ 

“  '  (n  + 


lV 


(uniaxial,  pure  tension) 


(uniaxial,  pure  bending) 


n=l 


\t  na  Z  /  'x  \ 

=  -  7  /  ~(n~~+~  ^  (J’t)  (equibiaxial ,  pure  bending) 

n-1 


Z  =  maximum  tensile  stress 
m 


B  =  Beta  function 


The  coefficients,  an,  in  the  preceding  series  are  determined  from 
material  specimens  by  fitting  a  least-squares  polynomial  to  the  fracture 
data : 


in? 

s 


N  +  1  - 
N  +  1 


i' m 


If  four-point  bending  tests  have  been  used,  for  example,  then  the  an 
coefficients  can  be  expressed  in  terms  of  the  bn  coefficients  of  the 
least-squares  fit  to  the  fracture  data. 

b 

a^  -  2(n  +  l)(2n  +  1)  — —  (uniaxial,  shear  insensitive) 


a 

n 


2(n  +  l)2 


n. 

V 


(uniaxial,  shear  sensitive) 


The  probability  of  survival  of  a  material  volume  element  in  a  triaxial 
stress  state  can  then  be  expressed. 


in? 

s 


r 


i 


n  T_  n-1 
na  CT  K  dk 
n  x  c  c 
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where  ax  is  the  maximum  principal  stress  (tension)  in  the  volume  element 
and  Q/4t r  is  determined  from  the  stress  state  that  exists  within  the 
element.  The  probability  of  dome  survival  is  then  determined  in  the 
same  way  as  described  for  Weibull’ s  fracture  theory. 


The  current  investigation  employed  a  simplified  version  of  the  above 
equations.  The  analysis  was  constrained  to  using  shear  sensitive  cracks 
and  the  critical  crack  strength  parameter,  M,  was  simplified  to 


where  m  and  aQ  are  the  Weibull  parameters.  The  parameters,  m  and  c,  are 
determined  from  the  equations  given  in  the  section  describing  the  Weibull 
fracture  theory.  The  parameter,  aQ,  is  determined  from  the  following 
relations : 


a 

o 


a 

o 


a 

o 


exp 


exp 


li 

I" 


£nV  -  C 


£n  -y  -  Zn(m  +  1) 


V 

In  —  +  2nm  +  £nB 

4 


(uniaxial,  pure  tension) 


(uniaxial,  pure  bending) 


(equibiaxial,  pure  bending) 


The  probability  of  survival  for  these  special  cases  then  becomes, 


>,nP 

s 


(uniaxial,  pure  tension) 


(uniaxial,  pure  bending) 


(equibiaxial,  pure  bending) 


and  the  probability  of  survival  for  a  volume  element  in  a  triaxial  stress 
state  is  found  by  numerical  integration  of  the  expression, 


4nr 


m(m  +  1) 


c 


1 


dK 

c 
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where  2/ 4rr  is  evaluated  for  the  stress  state  existing  within  the  volume 
element.  The  probability  of  dome  survival  is  found  by  the  method  pre¬ 
viously  described  in  the  section  on  the  Weibull  fracture  theory. 

The  simplified  method  described  above  reduces  to  Weibull's  method 
for  uniaxial  stress  states,  yet  retains  crack  orientation  effects  on  the 
probability  of  fracture  in  biaxial  and  triaxial  stress  states.  In 
addition,  the  least-squares  curve  fit  of  the  material  specimen  fracture 
data  is  identical  to  that  used  in  the  Weibull  theory  for  uniaxial  fracture 
tests.  Consequently,  the  values  for  the  parameters,  m  and  aG,  shown  in 
Table  2  for  uniaxial  fracture  data  for  MgF2  can  be  used  directly  in  the 
Batdorf  fracture  theory.  Application  of  the  simplified  method  to  equibi- 
axial  fracture  data  (see  footnote  7)  for  MgF2  leads  to  the  values  for  m 
and  aQ  shown  in  Table  3.  A  comparison  between  the  least-squares  fit 
approximation  and  the  equibiaxial  fracture  data  for  MgF2  is  shown  in 
Figure  18. 


TABLE  3.  Equibiaxial  MgF?  Fracture  Data. 


Temperature 
Volume,  a 

o 

Surface,  0 


7  5°F 

6,524  psi 
19,379  psi 
4.318 


SURFACE  DISTRIBUTED  CRACKS 

The  fracture  theory  developed  by  3atdorf  (footnote  5)  for  surface 
flaws  assumes  that  such  flaws  are  cracks  oriented  normal  to  the  surface. 
The  effects  of  crack  orientation  on  probability  of  fracture  in  a  biaxial 
stress  state  are  also  included  in  the  theory.  The  following  expression 
is  used  to  determine  the  probability  of  survival  for  a  small  surface 
area  element  in  a  uniform  biaxial  stress  state: 


-  -  iA/Md  (?) 
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FIGURE  18.  MgF2  Equibiaxial 
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where 


M  =  probability  that  a  crack  of  critical  strength  0C  (or  less)  is 
present  in  the  surface  area  element  (AA) 


—  =  probability  that  the  crack  is  oriented  in  a  direction  such 

that  the  stress  normal  to  the  crack  plane  exceeds  the  critical 
crack  strength 


a 

c 


=  critical  crack  strength 


The  probability  that  the  crack  is  oriented  in  a  direction  that  will 
result  in  fracture  is  determined  by  setting  *  -  0  in  the  expression  for 

cos^ 9  in  the  section  on  the  volume-distributed  crack  theory. 

c 


sin  6 

c 


(shear  insensitive  cracks) 


sin  9 

c 


(shear  sensitive  cracks) 


Since  yj  must  be  equal  to  twice  9 

c 


(Figure  19) , 


u> 

T 


1 


(shear  insensitive  cracks) 


('shear  sensitive  cracks') 


(equibiaxial  limiting  value) 


FIGURE  19.  Surface-Distributed 
Crack  Coordinate  System. 
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If  che  critical  crack  strength  parameter,  M,  is  again  expressed  in 
series  form, 


=  Y  a  an  =  a  an  Kn 
/  j  n  c  /  j  n  x  c 


n=l 


n=l 


then  the  probability  of  survival  can  be  determined  analytically  for  the 
following  cases: 

Shear  Insensitive  Cracks 


ZnP  = 
s 


ZnP  = 
s 


-  It  [H  +  B  I) 


n=l 


(uniaxial,  pure  tension) 


n=l 


ZnP 


-  A  >  a  0 

/  j  n  m 

n*=l 


(uniaxial,  pure  bending) 
(equibiaxial ,  pure  bending) 


*,nP  » 
s 


ZnP  - 
s 


?nP  - 
s 


Shear  Sensitive  Cracks 

(uniaxial,  pure  tension) 

\) 

(uniaxial,  pure  bending) 


n=l 


it  /  ^  n  +  1 


E 

n=l 


a  B 
n  m 


-  aS  ^  a  Jn 

/_  v  n  m 


(equibiaxial,  pure  bending) 


n=l 


The  an  coefficients  are  again  determined  from  a  least-squares  fit  to  the 
fracture  data.  However,  for  the  surface  distributed  crack  theory,  a 
solution  can  be  obtained  for  any  biaxial  stress  state*  from  the  following: 


ZnP 

s 


Shear  Insensitive  Cracks 


a 

n 


x 


I 

n 


*See  Appendix  B 
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where 


1  .  -1  /  -2cx  -  b 

I  =  —  sin  /  7  ■■■  — 

°  x  C  Wb  -  4ac 


I,  *  -  -r—  I 
1  2c  o 


(K  £x<  1) 


(3b2  -  4ac) 


3  \4c"  16c3/  ° 


j  =  _  (2n  -  1)  b_  l  _  (n  -  1)  a  t 

n  n  2  c  n-1  n  c  n-2 


a  =  -K 


b  -  1  +  K 


c  =  -1 


Shear  Sensitive  Crack 


;.nP  =  .  y  a  G 

s  /  i  n  x  n+1 


whe  re 


Gq  =  «k) 

Gi =  i 


G,  =  e(k) 


G3  *  4  t1  + 


:4  =  -  3  K‘<(k)  +  |[i  +  K“ ]£ (k) 
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G2n+4 


(jn  +  1) 

(2n  +  3)  y  2n 


C2n  +  2) 
(2n  +  3) 


[  1  +  K.*"  ]  G 
v 


2n+2 


G 


2n+5 


(n  +  1) 
(n  +  2) 


K2gi  ii 

v  2n+l 


(2n  +  3) 

(n  +  2) 


[1  4-  K2]G 

y 


2n+3 


o  9 

k“  =  1  - 

y 

K(k)  =  complete  elliptic  integral  of  the  first  kind 

e(k)  =  complete  elliptic  integral  of  the  second  kind 

If  the  above  equations  are  derived  in  terms  of  the  two  Weibull  para¬ 
meters,  m  and  0Q ,  for  shear  sensitive  cracks,  then  the  critical  crack 
strength  parameter,  M,  becomes 


The  parameters,  m  and  c,  are  determined  from  the  equations  describing 
fveibul!  fracture  theory,  and  the  parameter,  a  ,  is  decermined  from  the 
following  expressions: 


a  =  exp  -  -  C 


0  =  exp 

o  r 


(uniaxial,  pure  tension) 


=  exp 


xp  ! -■  f  In  /lW  4*  y- r  \  -  C  |  (uniaxial,  pure  bending) 

|  m  l  \  ™  +  1/  JJ 

h-t1'  i) 


nA  -  <,n 


(equibiaxial,  pure  bending) 


The  probability  of  survival  for  these  cases  then  becomes, 

m 


in?  -  -  A  f  — 
s 


(k) 


In?  =  -  |  LW  + 
s 


LH  /£_\™ 

b  +  iIW 


In? 


TA 


-V) 


(uniaxial,  pure  tension) 


(uniaxial,  pure  bending) 


(equibiaxial,  pure  bending) 
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and  the  probability  of  survival  for  a  surface  element  in  a  biaxial  stress 
state  is  found  by  numerical  integration  of  the  expression, 


where  co/tt  is  evaluated  for  the  stress  state  existing  wit!',  in  the  surface 
area  element.  The  probability  of  dome  survival  is  then  det erniu^u  ‘.:v 

the  method  previously  described. 

This  simplified  method  reduces  to  Weibull's  method  for  uniaxial 
stress  states  and  retains  the  effects  of  crack  orientation  in  biaxial 
stress  states.  Values  for  the  parameters,  m  and  z  ,  for  ecuibi axial 
fracture  data  are  shown  in  Table  3. 


RESULTS 


The  volume  and  surface  distributed  crack  theories  of  both  Veibull 
and  Batdorf  (simplified  version)  were  applied  to  the  thermal  stress 
distributions  in  the  full-scale,  MgFo  domes  that  were  subjected  to  the 
simulated  free  flight.  A  comparison  between  the  predicted  and  observed 
probability  of  dome  :r -enure  --s  a  function  nr  rime  offer  tie  s inu. 
launch  is  shown  in  Figures  20  through  2i.  All  of  the  domes  retain- 2 
their  geometric  shape  throughout  the  simulated  free  flight,  even  after 
dome  fracture.  Mote  that  two  curves  representing  predicted  probahilirv 
of  dome  fracture  appear  in  each  figure.  These  tw<  curves  represent 
probability  of  dome  fracture  for  the  maximum  x~.d  minimum  thermal  tr_s~ 
levels  that  were  observed  at  a  given  t  ima  in  tn^  simulate:  free-'.  :  f 
test  series. 


DISCUSSION 


Examination  of  the  fracture  statistics  that  are  available  for  MgFi 
shows  that  (1)  the  surface  scratches  on  the  uniaxial  flexural  test 
specimens  are  much  larger  than  on  a  dome,  and  <  the  tvmr^racuro 
dependence  of  the  flexural  strength  for  the  oc.jibiaxial  t^st  specimens 
is  not  known.  Consequently,  one  would  expect  chut  the  combination  of 
large  surface  scratches  and  the  effects  of  test  specimen  temper.-; tare 
would  result  in  higher  than  actual  nrobabi  1  i u it* s  of  done  fracture  for 
the  fracture  statistics  obtained  from  the  uniaxial  flexural  tests. 

Since  the  flexural  strength  of  'gF?  decreases  with  temperature,  use  of 
fracture  statistics  derived  from  eouioiaxial  flexural  tests  at  room 


* 
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temperature  will  tend  to  underpredict  the  probability  of  dome  fracture. 
The  optimum  fracture  theory  for  MgF?  domes  should  therefore  overpredict 
the  probability  of  dome  fracture  when  based  on  the  uniaxial  fracture 
data,  and  underpredict  the  probability  of  dome  fracture  when  based  on 
the  eauibiaxial  fracture  data.  It  is  evident  from  the  results  shown  in 
Figures  20  through  24,  that  only  the  surface  distributed  crack  theory  of 
Batdorf  satisfies  this  criterion. 


CONCLUSIONS 

1.  Fracture  statistics  obtained  from  flexural  tests  of  small  MgFo 
specimens  can  be  used  to  predict  the  probability  of  fracture  of  full- 
scale  domes  in  an  aerodynamic  heating  environment. 

2.  The  surface  distributed  flaw  theory  of  Batdorf  provides  an  accurate 
means  of  predicting  the  probability  of  dome  fracture  for  the  material, 
MgF  2. 
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FIGURE  22.  Comparison  of  Batdorf  Volume  '"heorv  and  Dome  Fracture 
Data  (Uniaxial  Fracture  Statistics). 
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FIGURE  24*  Comparison  of  Satdorf  Surface  Theory  and  Dome  Fractur 
Data  (Equibiaxial  Fracture  Statistics), 
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Appendix  A 

BATDORF  VOLUME  DISTRIBUTED  FLAW  THEORY 


=  principal  stress  in  the  x  direction 
0 =  principal  stress  in  the  y  direction 


O’  =  principal  stress  in  the  z  direction 

0^  =  stress  resulting  from  the  principal  stresse 

normal  to  the  plane  of  the  crack 

a  =  component  of  a„  in  the  a  -  a  plane 
n  N  y  z 

9  =  angle  between  the  0^  axis  and  the  crack  normal 

=  angle  between  the  O  axis  and  a 

y  n 


Then  for  shear  insensitive  cracks. 


3  =  a  cos  d>  +  0  sin  <|> 

n  v  z 


7 


=  a  [K  -  (K  -  K  )sintap] 

X  V  V  2 

2  0 

aXT  =  a  cos  9  +  a  sin”9 

N  x  n 

=  0  i  [K  -  (K  -  K  )sin“-' 

X  (  v  v  2 

!  +  [(1  -  Kv) 

+  (K  -  K  ) s in“ b ] cos“9 > 
v  z  J 

and  for  shear  sensitive  cracks 


cos  +  a  sin 


=  a  [K.  -  (K  -  K  ) sin  £] 

x^v  v  z 

2  2 

o.T  =  a  cos  9  +  O'  sin“9 
N  x  n 


7  i[K  -  (K  -  K  )  sin2>t  ]  +  (1  -  K 
x(v  v  z  v 


+  (Kv  -  Kz)  sin“vp]cos29| 


direction 
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TN  =  (ax  “  an)sin  e  cos  S 

=  a  [  (1  -  K  )  +  (K  -  K  )sin^  $]sin  9  cos 
x  v  v  2 


2  2  ,  2 
°s  =cn  +  tn 


where 


=  shear  in  the  crack  plane 

=  effective  stress  normal  to  the  crack  plane 
k  e 

K  -  a  /a 

y  y  x 

K  =  a  /a 

2  Z  X 


The  limiting  angle,  0C,  for  which  O ^  >_  Gc  or  >_  ac  is  found  by  setting 

rr  ~  rr  rr+ .  =  /r 


aN  “  acJ  or 


-  K  -  [K  -  (K  -  K  )  sin“1>  ] 

2  c  v  v  z 

cos  0  -  - -  -  — — 

1  -  [K  -  (K  -  K  )  sin“p] 

V  V  3? 


(shear  insensitive) 


~  K  -  [K  -  (K  -  K  )  sin“ }  1  “ 

2~  c  v  v  z 

cos  0C  =  - - — 


1  -  [K  ^  (K  -  K)sin; 
v  v  z 


(shear  sensitive) 


where 


K  =  0  /a 
c  c  x 


The  solid  angle  over  the  stress  ellipsoid  for  which  _>  Cc  or  C\j  uc 
is  then  determined  from  the  expression. 


=  8 


or 


4tt 


•o 

»\3 

C 

sin  0  d0  -  4tt  -  8 

Jo  - 

o  « 

T 

o’ 

r“ 

1  -  - 

1  cos  A  d6 

"T 

c 

COS  0  d0 

c 


o 
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The  following  special  cases  have  been  worked  out  for  both  the  shear 
insensitive  and  shear  sensitive  cases: 


Shear  Insensitive  Cracks 


n 

*  0) 

(uniaxial ) 

(K  =  1, 

y 

K  =  0) 
z 

(equibiaxial) 

Shear  Sensitive  Cracks 

%=  =  1  -  K 

c 

n 

>, 

u: 

Nrf' 

=  0) 

(uniaxial) 

SL  =Ji-K2 

4tt  f  c 

(K  =  1, 

V 

K  =  0) 
z 

(equibiaxial) 

Expressions  for  the  quantity,  Q/4tt,  have  currentlv  been  derived  onlv 
for  shear  insensitive  cracks  in  a  general  triaxial  stress  state.  The 
expressions  derived  for  shear  insensitive  cracks  are  as  follows: 

(K  >  K  >1) 
v  —  c  — 
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(o  <  k  <  :<  ) 

—  c  —  z 


4tt 


1 


where 


AQ(3,k)  is  Heuman's  Lambda  function 
Using  the  approximation, 

A  (3»lO  =  (1  -  k“)sin  p  +  —  k^S 

the  above  expressions  for  2/4^  become, 

(K  <  K  <  1) 
v  —  c  — 


(K  <  K  <  K  ) 
z  —  c  —  y 


- 

(i 

i 

u 

u 

1 

fl  -  K 

1  4. 

i  __  _ _ _ 

V  c  3 

t/  c 

47  = 

1  (i 

-  K  )  (K  -  K  ) 

X  1  -  K 

. 

C  V  z 

f  z 

2  (i  *  Kv)(Kc  '  V  .  -1 

TT  (1  -  K  )  (K  -  K  )  bln  \  1  -  K 
c  V  z  ^  z 


(0  <  K  <  K  ) 

—  C  —  2 


*  L 

readily  reduced  to  the  results  obtained  for  the  special 
cases  of  uniaxial  and  equibiaxial  stress  states.  The  corresnonding 
expressions  for  shear  sensitive  cracks  involve  rather  messy  hyperellipti 
integrals  and  are  probablv  best  treated  via  numerical  integration 
techniques . 


4tt 

which  can  be 
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form, 


If  che  critical  crack  strength  parameter,  NT,  is  expressed  in  the 


w  «  n  \  '  n  t,n 

N=>  a  a  *  >  a  a  K 
n  c  n  x  c 


n=l 


n=l 


Then  the  following  expressions  can  be  developed  from  the  special  case 
relations  for  Q/4tt: 


Shear  Insensitive  Cracks 


In?  =  -  V 
s 


n  m 


(2n  +  1) 


n*l 


(uniaxial,  pure  tension) 


2nP  = 
s 


V  \  1 _ n"m 

2  2^  (n  +  1)  ( 2n  =  1) 

n=  1 


(uniaxial ,  pure  bending.) 


cnP  * 
s 


-  y y  B  |  (equibiaxial  ,  pure  bending) 


n=l 


Shear  Sensitive  Cracks 


ZnP  = 
s 


n  =  J 


tnP  = 
s 


In? 


V 

? 


n«  1 


(n 

*  1) 

_n 

a 

n  m 

(n 

+  ir 

n 

(n 

+  1)  a 

(uniaxial ,  pure  tension) 


(uniaxial,  pure  bending) 


■  (n  +  TT  anam  B  (  2 ’2  )  <e(luibiaxial  .  Pure  bending) 


n=l 


where 

0  =  maximum  tensile  stress 

m 

3  =  Beta  function 
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The  coefficients,  an,  in  the  series  above  can  be  determined  from 
material  specimens  by  fitting  a  least-squares  polynomial  to  the  fracture 
data. 


>>nP 

s 


l 

n 


f  g  ±  1  -  nl 
L  N  +  1 


s 

S  —  1 


b.a1 
i  m 


where 

N  =  total  number  of  specimens  tested 

n  -  number  of  specimens  that  failed  at  a  maximum  tensile  stress, 

a  ,  or  lower 
m 

b.  =  least-squares  polvnomial  coefficients  for  fracture  data 

l 

If  the  four-point  bending  tests  have  been  used,  for  example,  then 
the  an  coefficients  can  be  expressed  in  terms  of  the  bn  coefficients  of 
the  least-squares  fit  to  the  fracture  data. 


b 

a  =  2(n  +  l)(2n  +  2)  ~ 
n  V 


(uniaxial,  shear  insensitive) 


a 

n 


2(n  +  1) 


.y  b 
2  n 


(uniaxial,  shear  sensitive) 


The  probability  c  survival  of  a  material  volume  element  is  then 
expressed  as 


for  triaxial  stress  states.  The  maximum  principal  stress,  7X,  can  be 
determined  from  a  finite-element  stress  analysis  of  the  geometrv  in 
question  and  the  an  coefficients  are  determined  from  fracture  test  data 
obtained  in  a  more  convenient:  stress  state  (e.g.,  uniaxial1).  The  proba¬ 
bility  of  survival  for  the  geometry  in  question  is  then  equal  to  the 
product  of  the  survival  nrobabil i t ies  for  the  individual  \  1  ume  elements 
in  the  finite-element  model. 
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Appendix  B 


BATDORF  SURFACE  DISTRIBUTED  FLAW 

If  8  =  constant  =  0,  in  the  expression  for  • 


The  probability  of  survival  of  an  element  of  suri 


THEORY 


(Appendix  A) ,  then 


(shear  insensitive) 


(shear  sensitive) 


(shear  insensitive) 


(shear  sensitive) 


ace  area  is  c;iven  by 


(shear  insensitive) 


(shear  sensitive) 
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If  che  critical  strength  parameter,  M,  is  expressed  as  the  series 


The  integral  for  the  shear  insensitive  case  is  of  che  form. 


where 

v  s  a  +  bx  +  cx” 

and 

a  =  b  =  (1  +  ) ,  c  =  -i 

Evaluating  these  integrals. 


N'WC  TP  t>22o 


and 


I  =  -  (2n_- JL>.  b_  1  _  (n_-_U  a  z 

n  n  2  c  n-1  n  c  n-2 


Therefore, 


-  nP 


-TE 

n=l 


a  a  I  dA 
ax  n 


(shear  insensitive) 


The  corresponding  integral  for  the  shear  sensitive  case  reduces  to 
the  form 


G  , ,  =  G  = 
n+1  m 


(m-1) 


r  2  2 1  2 

I  1  -  k  sin  :<  I  dx 


where 


x  =  sm 


■i is, 


9  2 
k~  -  1  -  K 

y 


t  =  K" 


Evaluating  these  integrals 

co  -  <(K) 

Gi  -f 

G,  =  (k) 

G  »  r  (1  +  K2) 


and 


G,  =  -  1  K“V(k)  +|(1+  K2)e(k) 
4  3  y  3  v 


2n-K 


(2n  +  1)  2  .  (2n  +  2)  ,  2 

(2n  +  3)  vy  *2n  (2n  +  3)  y  2n+2 


,  _  (JL±J1  K-G  x  1  (1  +  /- 

2n+5  (n  +  1)  y  2n+l  2  (n  +  2)  v  y  2n+* 


Therefore 


-"nP 


IE 

11=1 


a  anG  dA 
ti  x  n4-! 
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